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Abstract 

Since Kolmogorov proposed his phenomenological theory of hydrodynamic turbulence in 1941, 
the description of mechanism leading to the energy cascade and anomalous scaling remains an 
open problem in fluid mechanics. Soon after, in 1949 Onsager noticed that the scaling properties 
in inertial range imply non-differentiability of the velocity field in the limit of vanishing viscosity. 
This observation suggests that the turbulence mechanism may be related to a finite-time singularity 
(blowup) of incompressible Euler equations. However, the existence of such blowup is still an open 
problem too. In this paper, we show that the blowup indeed represents the driving mechanism 
of inertial range for a simplified (shell) model of turbulence. Here, blowups generate coherent 
structures (instantons), which travel through the inertial range in finite time and are described by 
universal self-similar statistics. The anomaly (deviation of scaling exponents of velocity moments 
from the Kolmogorov theory) is related analytically to the process of instanton creation using the 
large deviation principle. The results are confirmed by numerical simulations. 
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I. INTRODUCTION 



Describing the mechanism of developed turbulence for the 3D Navier-Stokes equations 
remains an important open problem in fluid mechanics. It encompasses various questions, 
and in this work we address the anomalous statistics of velocity moments in inertial range and 
the dissipation anomaly (existence of finite dissipation in the inviscid limit). These questions 
remain the hot research topic since Kolmogorov presented the phenomenological theory of 
inertial range in 1941 Jj|. This theory of isotropic homogeneous turbulence leads to the 
power-law dependence of velocity moments on spatial scales, providing the scaling exponents 
Cp = p/3 obtained on dimensional grounds. The exact scaling exponents deviate from the 
Kolmogorov theory. These deviations, called the anomalous corrections, are universal and 
become large with increasing p. Though a lot of knowledge is available now on the described 
anomalous phenomena, their mechanism is still not well understood j^, 3|. 

In 1949, Onsager 4| related scaling properties of turbulent flow in inertial range with the 
regularity of solutions obtained in the limit of vanishing viscosity. He conjectured that the 
anomalous turbulent dissipation requires the limiting velocity field to be non-differentiable 
with the Holder continuity exponent h < 1/3. This conjecture was proved later 0, 6]. 
Irregularity of inviscid solutions allows considering the flow as a multifractal set with a con- 



tinuous infinity of dimensions [2|, tJ, which explains the nonlinear shape of scaling exponents 
( p . Development of the theory of turbulence in this way has the fundamental obstacle. It is 
the problem of blowup, i.e., the formation of a finite-time singularity in the incompressible 
3D Euler equations from smooth initial data of finite energy. So far, the existence of blowup 
remains an open problem 8|. 

Simplified models help in understanding the turbulence phenomena. In this respect, the 
Gledzer-Ohkitani-Yamada (GOY) shell model of turbulence was successful in describ- 

ing several nontrivial properties including the inertial range with anomalous dissipation and 
scaling. Shell models represent the dynamics in terms of characteristic (shell) velocities 
corresponding to a discrete set of wavenumbers increasing in geometric progression, and 
allow reliable numerical simulation at very high Reynolds numbers. The Sabra shell model 



proposed in Uj is characterized by improved regularity in the inertial range. Despite of 
large effort v\ . the theory of turbulence for shell models, which would follow directly from 
the model equations and describe the observed statistics, is not yet accessible. On the other 



blowup in the Sabra shell model has self-similar universal structure 
of such structure to the statistics of turbulence was discussed in 
turbulence 
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hand, the problem of blowup was recently formalized 13] and understood fl^ . The 



16 1 . Possible relation 

V 

181 . Cascade models of 



191 ] represent the extended version of shell models, where each shell is described 



by a large 
dynamics 



though fixed) number of variables. These models lead to anomalous intermittent 



20|, and the universal sel 



cid cascade model 



2l|. See also 
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similar blowup was observed numerically in the invis- 



22M24I| for other numerical observations of self-similar 



blowup in shell models. 



In this paper, we establish a direct link between the blowup and the turbulent dynamics in 
inertial range for the Sabra shell model. We show that the blowup-like structures dominate 
the turbulent fluctuations and can be described as a "gas" of instantons. The instantons 
(coherent structures of shell velocities, which traverse the inertial range in direction of large 
wavenumbers) are represented and analyzed in terms of velocity local maxima. The striking 
property of instantons is that they maintain the universality and self-similarity of blowup, 
though with slightly different scaling exponents and in statistical sense. This statistical 
universality of instantons was observed earlier in 25[ . Then we show that instanton creation 
is the main process driving the inertial range dynamics. This allows deriving the probability 
density function (PDF) for instanton amplitudes explicitly in terms of the anomalous scaling 
exponents ( p by using the large deviation principle. The obtained results fully agree with 
numerical simulations and are also confirmed analytically for a class of instanton creation 
models. Finally, we discuss some qualitative changes in the turbulent regime, which occur 
with a change of model parameter. 



The paper is organized as follows. Section [XX] introduces the Sabra shell model. The 
blowup universal properties in the inviscid model are described in Section [TTTJ In Section HVl 
we consider statistics of maxima of velocity amplitudes and introduce a way to identify 
the instantons. Section |V] describes the universal self-similar statistics of instantons. In 
Section IVI1 we find universal expressions for PDFs of instantons using the large deviation 
principle. Section IV HI presents the analytical theory for a specific class of instanton creation 
models. Section [Villi describes a different turbulent regime, which is dominated by a single 
blowup. The results are summarized in Section [DO 
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II. MODEL 



In shell models of turbulence, the Fourier space is represented by a series of shells n = 
0, 1, 2, . . . corresponding to wavenumbers k n = X n with A = 2. We consider the Sabra shell 
model 11] 



= i[k n+1 u n+2 u* n+1 - (1 + c)fc n ii n+ iu*_i 
at (1 ) 

-ck n -iu n -iu n -2] - vk 2 n u n + /„, 
where u n is the complex shell velocity, which can be understood as the Fourier component of 
the velocity field at the shell wavenumber k n , v > is the viscosity, and c is the parameter 
controlling nonlinear coupling of the shells. The terms f n model external forces at large 
scales and, thus, they are usually restricted to the first few shells. The inviscid system 
with no forcing (y — f n = 0) conserves the energy E = ^Yl n \ u n\ 2 - The second quadratic 
invariant H = 'Yl in c~ n \u n \ 2 is associated with the helicity for c = —0.5 when c~ n = (—l) n k n . 
Additionally, there are four symmetry transformations 

t^t-to; (2) 

u n ^ e ie "u n , 6 n = + # n _ 2 ; (3) 
t\-tt/a, u n ^au n ] (4) 

U n XUn+l- (5) 



Here Eqs. (J2]) and ([3]) can be associated with the time and physica 



Eqs. OH) and <^ correspond to the time and space scaling, see [12]. 



space translations, while 



III. BLOWUP IN INVISCID MODEL 



Let us consider solutions u n (t) with the finite norm |n|i < oo defined as 

t \ 1/2 

i=(E^H 2 • ( fi ) 

v n J 

The blowup represents a singularity given by 

— > oo as t — > t~, (7) 
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which develops in finite time t c < oo from initial condition of finite norm 131 ] . Note that 
the singularity is described by the norm, while each particular shell speed u n {t) remains 
finite and smooth. This reflects the fact that the shell model corresponds to dynamics in 
the Fourier space, where the condition like (|7|) implies the divergence of velocity derivatives 
in physical space, i.e., infinite vorticity. The blowup is only possible in the inviscid shell 
model, and the uniqueness of solution is insured only for t < t c 13]. 

Let us consider the inviscid model with vanishing forcing terms, v — f n = 0. Then we 
write Eq. ([T]) as 



dt 



^ = N n [u'}, u' n = ik n u n , 



(8) 



with the quadratic nonlinearity 



N n [u'} = -A- 2 < +2 < +1 + (1 + cK +1 <-i " c\ 2 u' n ^u' 



n-2- 



(9) 



Following the approach suggested by Dombre and Gilson [14( (see also [3]), we consider the 
renormalized time r and shell speeds w m introduced as 



t = t + / exp 
'o 



exp 



A(r')dT' 



A(r")dr'' 



dr' 



(10) 



where r = corresponds to the initial time t , and A{r) is specified below. It is straightfor- 
ward to check that 

rim 

N n [w]-Aw n , (11) 



dr 

where iV n [w] has the form written in terms of w n instead of u' n . One can also check that 
Eq. (1TTT) conserves the sum ^2 \ w n \ 2 if we choose 



A(r)= Re 



W r , 



(12) 



The idea of the above transformation is that Eq. f TTTj) admits an asymptotic traveling 
wave solution of the form [l^] 

w n (r) =W(n-sr), (13) 

where s is the wave speed and W(!;) is a function vanishing as £ — > ±oo. This traveling wave 
exists for a large range of shell model parameter c, and it is determined up to symmetries 



induced by Eqs. For 
(fT0|), solution (H3D yields Q, 



u n (t) = -ik-y°U(k z n °(t-t c )), 



where 



;he original shell speeds u n {t) related to w n (r) by Eqs. ([8]) and 

(14) 

(15) 
(16) 
(17) 



U(t-t c 



1 



exp 

l/s 



log A J 



A{r')dr' W(-st), 
A(r)dr, y =l- z , 
dr'. 



POO 




/ exp 


- [ T A(r")dr" 


/o 


Jo 



t c — to 

If yo < 1, then Eq. (1141) describes the asymptotic form of blowup at finite time t c < oo. In 
this asymptotic form, y is the universal scaling exponent independent of initial conditions, 
and the function U (t) describes the universal self-similar shape of the blowup given up to 
the scaling symmetry of the Sabra model. The equality yo + zq = 1 reflects the dimensional 
relation t n — t c oc (v n k n )~ l , where v n = max ( |w n (^)| an d t n is the corresponding time. For 
details of the derivations and the rigorous theory, which associates the traveling wave (Tl3l) 
with a fixed-point attractor of the Poincare map, see [l^j. 

As an example, let us consider the case c = —0.5. Solution w n (r) of the renormalized 
system (TTTT) for real initial conditions is shown in Fig. [Th. One can clearly see the formation 
of traveling wave solution (fT3l) . Solution for the original shell speeds u n (t) is presented in 
Fig. [lb, which blows up at finite time t — > t c given by Eq. (fT7j) . Using Eqs. (jT5"l) and (fl6l) . 
we compute the scaling exponent yo = 0.281 and the function U(t). The bold green curves 
in Fig. [Tb show the asymptotic self-similar solution ( T141 for the blowup, and one can readily 
see the convergence. Numerical analysis confirms asymptotic stability of the traveling wave 
solution in Fig. due to both real and complex perturbations. As we already mentioned, 
this implies that Eq. (TMj) provides the universal asymptotic form of blowup. 

Similar traveling wave solutions exist for c < —0.092. The corresponding scaling exponent 
yo and function U(t) are shown in Fig. [2j The function U(t) is monotonous for c < —0.671, 
possesses a single extremum (maximum) for —0.671 < c < —0.139, and has several extrema 
for —0.139 < c < —0.092. At c = —0.139, we have U(0) =0 and the scaling exponent 
attains the minimum y Q = 0. This fact can be understood using the energy conservation 
argument. Indeed, y in Eq. fTT4|) cannot be negative, otherwise the shell speeds and the 
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FIG. 1. (Color online) (a) Traveling wave formation in the dynamics of renormalized inviscid Sabra 
model. Shown are the curves w n (j) with n = 0, 1, . . . increasing from the left to the right, (b) The 
corresponding dynamics of shell speeds u n (t) for n = 2, 3, . . .. Bold green (light gray) curves show 
the universal self-similar asymptotic form of blowup. Red squares indicate the correlated sequence 
of maxima v n = maxt |u n (t)|. 

energy would grow infinitely. In the case x/q = all the energy is transported to large shells 
as t — » t~ , so that no energy remains in each shell at the time of blowup, i.e., U(0) = 0. 

The real traveling wave solution ffTB"]) becomes unstable with respect to complex pertur- 
bations at the critical value c = —0.092. For c > —0.092 analysis of the blowup requires 



more sophisticated techniques, see [151 ] . which is beyond the scope of this paper. 



IV. INSTANTONS IN INERTIAL RANGE OF TURBULENT REGIME 

It is known that, for the parameter c = —0.5, the Sabra model with small viscosity 
(large Reynolds number) demonstrates chaotic intermittent behavior. Statistical properties 
of this system have much in common with the developed turbulence of the 3D Navier-Stokes 



equations [11] . In particular, it possesses a wide (increasing with the Reynolds number) 
inertial range of wavenumbers k n separating the scales influenced by forcing (small k n ) and 
the scales dominated by viscosity (large k n ). This inertial range is responsible for the energy 
cascade, i.e., to the flux of energy produced in the forcing range by external forces to the 
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FIG. 2. (Color online) (a) Dependence of the blowup scaling exponent yo on the Sabra model 
parameter c. Squares show the scaling exponents y of instantons. (b) The universal function U(t) 
in the asymptotic expression (fbi|) for c = —0.1 (upper curve, blue), c = —0.139 (dotted curve), 
c = —0.5 (solid black curve) and c = — 1 (lower curve, red). 



viscous range, where it is dissipated due to viscosity. Existence of a positive limit of mean 
dissipation rate for infinite Reynolds numbers constitutes the famous dissipation anomaly 
of turbulent hydrodynamic flows. 

The important quantitative characteristic of inertial range is given by the structure func- 
tions (velocity moments). In the inertial range, these functions depend on k n as power 
laws 

s p (k n ) = (KH ex fc-fc. (is) 

In this expression p is an arbitrary real number; traditionally, the computations are carried 
out for positive integer values of p. The scaling exponents ( p are universal, i.e., they are 
independent both of the forcing and viscosity. Figure E] presents the functions S p (k n ) in 
logarithmic coordinates for the Sabra model with c = —0.5. These results are based on 
direct numerical simulation of Eq. (pQ) with 40 shells, viscosity v = 10~ 14 and the constant 
forcing at the first two shells, /o = 1 + i and j\ = fo/2. One can clearly distinguish the 
forcing range corresponding roughly to the shells n < 5, the viscous range of shells n > 32, 
and the linear part in between indicating the inertial range. 

The phenomenological theory developed by Kolmogorov (K41 jl, 2|) predicts the linear 
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FIG. 3. (Color online) Solid black lines present the velocity moments S p (k n ) = (\u n \ p ) for p = 
1,...,6. Black circles determine the functions S'Jk n ) from Eq. (flU]) ; for better comparison of 
slopes, the graphs are shifted in vertical direction and only even n are shown. Red crosses show 
similar functions S' (k n ) computed for the local maxima v n corresponding to stable instantons only. 
Three types of structure functions determine equal slopes in the inertial range given by the scaling 
exponents — ( p . 

dependence ( p = p/3 for the scaling exponents. However, the exact scaling exponents ( p 
depend nonlinearly on p. This deviation from the K41 theory is called the anomaly. The 
scaling exponents are presented in Fig. HI The two exact values of scaling exponents are 
known. The first one is Co = since \u n \° = 1. The second exact exponent is (3 = 1, which 



is a necessary condition for the dissipation anomaly, see, e.g., U . The scahng exponents 
of the 3D Navier-Stokes turbulence are close to the ones given by the Sabra model [2|, IX XI ] . 

In this section we establish a link between the anomalous turbulent statistics and the 
blowup phenomenon for the Sabra shell model. The blowup analysis of the inviscid model is 
relevant in the inertial range, where viscosity is insignificant. However, there is an essential 
difference related to initial conditions. For the blowup considered in Section HH] finiteness of 
the norm requires decay of initial shell speeds faster than k~ l . This condition is violated 
in the inertial range of developed turbulence, which is characterized by the power-law decay 
(j!8p with d 0.39. We will see that this difference leads to the transformation of the 
blowup with universal self-similar asymptotic form to coherent structures with universal 
self-similar statistics. 
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FIG. 4. (Color online) Anomalous scaling exponents £ p computed for the velocity moments S p 
(black line) and for the functions S' (circles and crosses correspond to the sums over all maxima 
and over maxima from stable instantons, respectively). The red (gray) line Q p = p/3 corresponds 
to the phenomenological K41 theory. The dotted line shows the upper bound (|28p based on the 
instanton scaling. 

Identification of these coherent structures in turbulent regime is strongly facilitated, if 
we consider local maxima v n = max; of shell speed amplitudes. An extra subscript 

is necessary to index all the local maxima in shell n, but we will drop it for the sake of 
simplicity of notations. The new "structure" functions are defined as 

where the sum is taken over all local maxima v n observed for the speed amplitude 
during a large time interval < t < T. By a simple dimensional consideration, one finds 
that each local maximum v n has the characteristic time At n ~ (fc n f n ) _1 determining the 
time interval, where |u n (^)| ~ v n- F° r the velocity moment (\u n \ p ) = T~ x J Q T \u n \ p dt, this 
yields the contribution of order 

T- l vP n At n = (Tk^vfr 1 , (20) 

leading naturally to Eq. ( II 9p . Hence, the functions S' are expected to scale in the same way 
as S p in the inert ial range, i.e., 

S' p (k n ) oc ktf> (21) 
10 
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FIG. 5. (Color online) Typical dynamics of speed amplitudes |tin(^)| shown for the shells n = 
7, ...,24. Red squares mark correlated sequences of local maxima (instantons), which have the 
structure similar to the blowup in Fig. [TJ Shown are the instantons created in shells no = 7, . . . , 14. 



with the same scaling exponents as in Eq. ( I18p . This hypothesis perfectly agrees with the 
numerical simulations as shown in Figs. [3] and HJ 

The blowup in the inviscid shell model can be identified as the correlated sequence of 
maxima, which follow in increasing order of n and t, see Fig. [TJ Analogous correlated 
structures (called the instantons) are observed in the inertial range of shell models 
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26[, see Fig. Following [25|, we identify the instanton as a sequence of local maxima 
v n = max.t\u n (t)\ at times t n following in increasing order t no < t no +i < ••• < t ni . In 
this definition, no maxima of |u n (^)| or |u n+ i(t)| are allowed in the interval t n < t < t n+ \. 
Each instanton is created at some shell number uq and either reaches the viscous range or 
annihilates at a shell number rii in the inertial range. Using this rule, we group all maxima 
of velocity amplitudes into instantons, Fig. [5j 

As we already mentioned, an instanton can be viewed as a blowup deformed by the inertial 
range environment, in which it propagates. One can see from Fig. [5] that this deformation 
is caused, mostly, by interaction with adjacent instantons. Let N a u be the number of all 
maxima v n in a given shell n. Figure provides numerical values for the relative number 
N/N a u, where is the number of maxima in shell n corresponding to a specified type of 
instantons. Most of the maxima correspond to stable instantons, which reach the viscous 
range, i.e., in our simulation rii > 32. These instantons cover from 60 to 90% of the total 
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FIG. 6. (Color online) (a) The number N of selected local maxima in the shell n relative to their 
total number N a u. The bold black line corresponds to the maxima chosen from stable instantons. 
The thin red (lower) line corresponds to uncorrelated maxima (instantons of length 1 or 2). (b) 
Power-law scaling for the number M n of stable instantons created in shell n per unit time. The 
slope 1 — £i is shown by the dotted line. 

number of maxima in a given shell n (bold black line in Fig. Ek). Majority of the remaining 
maxima (about 20%) belong to very short instantons with n\ ~ uq, which can be considered 
as uncorrelated fluctuations. The instantons annihilating after traversing more than 2 shells 
but before the viscous range are rare. 

We see that the turbulent dynamics in inertial range of the Sabra model has the highly 
correlated structure, where the blowup plays a role of the driving mechanism. Another 
evidence supporting our observation is obtained if we compute the scaling exponents ( p for 
the functions (IT9]) . where only the maxima from stable instantons are included in the sum. 
These results are shown by crosses in Figs. [3] and HI The same scaling exponents as for the 
velocity moments ([I5|l are obtained (a tiny difference in Fig. H] is the same for all p and, 
thus, corresponds to a small change of the total number of maxima included in the sum). 
In the following analysis we will consider only the maxima belonging to stable instantons in 
the sum f fT9j) . 

Description of the inertial range in terms of instantons provides a new interpretations of 
the first scaling exponent (i. Let M n be the average number of stable instantons created in 
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shell n per unit time. Using Eqs. (119]) and ( )2T|) we have 

1 1 - 

S[(k n ) = ^l= ^M m oc (22) 

m=U 

where the first sum counts the maxima of stable instantons in shell n. It is easy to check 
that Eq. (|22)) implies 

M n <xk}r Cl . (23) 

We found that d determines the power-law scaling for the number of instantons created in 
shell n. For the Sabra model with c = —0.5, we have 1 — Ci = 0-61 in very good agreement 
with numerical data, Fig. [6b. 

The scaling exponent Co = is a simple consequence of the equality (|w n | ) = 1- However, 
this exponent gets nontrivial interpretation in terms of velocity maxima in Eq. (fT9l) written 
as 

where, as we showed earlier, At n ~ (A^n) -1 is the characteristic time associated with the 
maximum v n . Relation S' oc = 1 implies that the total fraction of time occupied by 
these maxima is finite for each shell, i.e., the stable instantons are dense in space-time. 

Numerical simulations for the model with the parameter values c = —0.4 and —0.6 were 
also carried out. The results are very similar to those presented in Figs. EMU which confirm 
our conclusion about the role of instantons as a principal elements of turbulent dynamics in 
the inertial range of the Sabra model. 

V. SELF-SIMILAR STATISTICS OF INSTANTONS 

The universal self-similarity of blowup Q14p is destroyed in the turbulent regime due to 
chaotic emergence and interaction of instantons, Fig. |5j The most striking property of the 
instantons is that they restore the blowup self-similarity in statistical sense in the inertial 
range. To observe this property, let us consider the functions 

4 no) (M = ^5>- n - n °' (25) 

(no) 

where the sum is taken over the local maxima belonging to stable instantons created in fixed 
shell uq. These functions can be viewed as effective velocity moments for the instantons born 
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FIG. 7. (Color online) (a) The functions R^ l °\k n ) in logarithmic coordinates demonstrating power- 
law scaling of instantons. Curves of the same color correspond to the instantons created in shells 
no = 13, . . . , 23. Different colors indicate different values of p = 1, . . . , 6 from top to bottom, (b) 
Graphs of the left figure collapse onto a single straight line when divided by p. The slope — y is 
shown by the dotted line. 

in a specific shell, and their graphs obtained numerically are shown in Fig. [7^ in logarithmic 
coordinates. One can clearly see that the functions obey the power-law scaling with 

exponents (slopes) independent of the initial shell number no- 

The next observation is that the slopes of the graphs in Fig. [7^ are proportional to p. 
This is shown in Fig. [TJd, where the functions (1/p) log 2 are plotted versus a number 
of shells n — n traversed by the instanton. All curves (after the vertical shift) collapse onto 
a single straight line of slope — y with y « 0.22. This implies the relation 

R^\k n ) = c C«o) A - ro An K fc -w An = n - n > 0, (26) 

with the universal value of scaling exponent y in the inertial range. The scaling exponent 
y 0.22 is different but close to the scaling exponent yo ~ 0.281 of the blowup, see Fig. 

The scaling rule in Eq. ( 1261) suggests the universal self-similarity of instanton statis- 
tics. Let us consider the probability density functions (PDFs) determining the probability 
Pn ,n(v)dv to sample a local maximum v = max t |w n (£)| belonging to the instanton created 
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FIG. 8. (Color online) Renormalized PDFs P no (v) of instantons found numerically for no = 17 and 
n = 17, . . . , 29 (thin black curves) and for no = 20 and n = 20, . . . , 29 (thin blue curves). Collapse 
of the graphs with fixed no onto a single curve confirms the self-similarity of PDFs in inertial range. 
Bold dotted curves show the PDFs determined by the large deviation principle. 

in shell no- The self-similarity for PDFs implies that the renormalized function 

PM = X^P^X^v) (27) 

does not depend on n in the inertial range. This hypothesis fully agrees with the numerical 
results as one can see in Fig. [Bk,, where the functions (127|) for different n collapse onto a 
single curve for fixed no = 17 or 20. The functions P nQ { v ) for different n are related by the 
large deviation principle, as we will show in the next section. 

We conclude that the instantons created in a given shell possess self-similar statistics. 
These instantons can be viewed as the blowup phenomena, which propagate to the viscous 
range interacting with each other. Interaction is an important factor which leads to a small 
but finite difference between the scaling exponent of the instanton y ps 0.22 and the scaling 
exponent of the blowup yo ~ 0.28. Similar results are obtained for the Sabra model with the 
parameters c = —0.4 and —0.6. The corresponding values of scaling exponents y are shown 
in Fig. EK. 
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VI. LARGE DEVIATION PRINCIPLE FOR INSTANTON DISTRIBUTIONS 



According to Eqs. 019p . (j25j) and (|26|) . an average contribution of a single instanton to 
the function S'(k n ) is proportional to kn 1 ^ ■ This yields an upper bound for the scaling 
exponents in Eq. ( 121]) as 

t p <l+(p-l)y. (28) 

The dotted line in Fig. H] represents the right-hand side of Eq. (128]) . Since the graph of ( p is a 
concave function j^J , we conclude that the instanton scaling exponent y does not determine 
any part of the ( p graph. In particular, y « 0.22 is larger than the slope of the ( p graph for 
large p (the numerical data provides the slope d( p /dp decreasing below 0.19). Therefore, the 
instanton scaling does not determine the scaling of high-order velocity moments, as it was 



181 ] (however, this becomes true for different values of the model parameter 



conjectured in 
c, as we show in Section IVHII) . 

The anomalous exponents C p arise in the process of instanton creation. In order to see 
this, we use relations (|T9|) . (1251) (1261) and find 

n 

S'p(kn) = k n x R p °\{k n ) 

no=0 
n 

no=0 

Then the coefficients are expressed from (1291) as 



(29) 



c 



(n) 
p-1 



k n S'(k n ) - X~^ 1)y k n _ 1 S'(k n _ 1 ). (30) 



In the inertial range, where the power-law scaling (I2ip holds, we have 

c p % « k^. (31) 

This relation was also confirmed numerically. We see that, due to the self-similar structure 
of instantons, anomalous scaling is attributed exclusively to the coefficients Cp™ ^ describing 
amplitudes of instantons created in shell uq. This property relates the inertial range anomaly 
with the process of instanton creation. 

Relation (|3~T1) allows finding the universal form of PDFs P n (v) in Eq. (T27|) as follows. 
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Using Eqs. (125]) -(ETt we obtain 

C p-1 = Rp-l{kn) = —'y]Vn 

(n) (32) 



COD 

M, 



roo poo 
n / V^P^dv = M n \ V^ X P n {v)dt 

Jo Jo 



where M n is the number of instantons created in shell n per unit time. Introducing the new 
variable a and function p(a) as 

1 v 

a = - log A — , p n {a) = p*nM n P n {v), (33) 
n 

where v* and p* are constant coefficients, we write Eq. (1321) in the form 

& = — log A / X npa Pn (a)da. (34) 
P* J 

Using Eq. fl3T|) . we find the power-law scaling rule for the integral in the right-hand side as 

J X npa p n (a)da oc klr Cp = X n{1 ^\ (35) 

In Eq. (|35l) the scaling exponent 1 — ( p is a smooth convex function of j) 6 R, Fig. HI 

27H29| to Eq. fl35|) . which states that p n {a) 



Hence, we can apply the Gartner-Ellis theorem 
has the asymptotic form 

p n (a) oc \~ nJ{a) = fc" J(a) (36) 

for large n, where the rate function J(a) is the Legendre transform of the function 1 — £ p , 
i.e., 

J(a)=pa-(1-Q, a = -^. (37) 

Expression ( 13"B"|) is called the large deviation principle. Note that the Gartner-Ellis theorem 
is formulated for A = e but one can easily check its validity for any A > 1. 

We verify Eq. (|36|) in Fig. HI where the black curves show the functions — (1/n) log A p n (a) 
found numerically using Eqs. (|27p and (133]) for n = 17, . . . , 26 and no = n. As expected, 
these graphs collapse onto a single curve given by the rate function J (a). The rate function 
represented by the red dotted line was computed using the Legendre transform f[3T|) for 
the scaling exponent ( p in the interval —2 < p < 10. In numerical computations, it was 
important to choose good values of the constants and in Eq. ( 133]) in order to achieve 
better convergence. 
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FIG. 9. (Color online) The functions — (1/n) log^ p n computed numerically for n = 17, . . . , 26 (thin 
black lines) are compared with the rate function J (a) (dotted red line). 

The final result of our derivation is obtained by substituting Eq. fl36l) into fl33|) as 



where we dropped the factor n representing a logarithmic correction for the first expression. 
Note that the asymptotic form given by the Gartner-Ellis theorem in Eq. ( 1361) is understood 
as n -1 log A p n (a) — > — J (a) in the limit n — > oo. Recall that the same limit of large n is used 
in the definition of inertial range, which corresponds to shell numbers far from the forcing 
range. Thus, Eq. fl38|) is valid in the inertial range. This statement is confirmed numerically 
in Fig. |H1 where the asymptotic PDFs given by Eq. ( )38|) are shown by the dotted red curves 
for the shells 17 and 20 (with the constant factors properly adjusted). 

We showed that the PDFs of instantons in the inertial range have the universal self-similar 
form ( 1381) related to the anomalous scaling exponents by Eqs. ( 1371) . Thus, the instantons 
satisfy the large deviation principle leading to the inertial range anomaly. The presented 
analysis has much in common with the phenomenological model of multifractality j^jj. In 
this model, it is assumed that the velocity field can be decomposed into fractal subsets 
with different scaling properties, and the fractal dimensions are related to the anomalous 
exponents ( p by the Legendre transform. However, the fractal subsets in the multifractal 
model are hard to define and detect numerically or experimentally, as well as to justify 
their appearance. On the contrary, the presented approach based on the study of instantons 
is related to the analytical theory of blowup and is supported by the detailed numerical 




(38) 



18 



analysis. 



VII. INSTANTON CREATION MODEL 

In this section, we introduce a phenomenological model for instanton creation, where the 
large deviation principle can be derived analytically. As one can see in Fig. \5\ an instanton 
traveling through the inertial range leaves a trace (energy) in all the shells it passed through. 
Due to the asymptotic stability of blowup mentioned in Section IHIt this energy "feeds" a 
series of newly created instantons in different shells. This process leads to formation of a 
"gas" of instantons, which is dense in space-time and carries the energy from the forcing to 
the viscous range. As the viscosity plays no role in this process, the dissipation anomaly 
becomes a natural consequence of the described behavior. 

In this phenomenological picture, instantons create other instantons. A simple statistical 
model of the creation process can be developed as follows. We assume that an instanton, 
which reaches the shell n with the amplitude v n = 1, creates in average (p(v)dv new in- 
stantons of amplitude v n = v in this shell. Here (p(v) is the creation rate function, which 
is assumed to be universal, i.e., independent of n. For an instanton of arbitrary amplitude 
v n = v', the density of created instantons is given by 

MS)*- (39) 

as it follows from the scaling symmetry of the Sabra model. As before, we consider only 
stable instantons, which cover up to 90% of all shell oscillations (Fig. [6^), and disregard 
other types of fluctuations. 

Using the definitions of Section |V] the distribution of instanton amplitudes is described 
by the product 

M n P n , n (v)dv, (40) 

determining a number of instantons with maxima v n = v created in shell n per unit time. 
Here P n>n is the PDF of such instantons and M n is the total instanton creation rate in shell 
n. Distribution of maxima v n = v corresponding to the instantons created in previous shells 
n < n is found similarly as 

71-1 

M no P n0tn {v)dv. (41) 

"0=0 
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Using expressions ( 140]) and (j41j) . the instanton creation principle described by Eq. (|39|) yields 



M n P n , n (v) = 
Using Eq. fl27|) we write this expression as 



n-l 
_n =0 



n ,nV 



v' \v' / 



Let us compute the quantity (1521 represented as 



v p M n P n (v)dv. 



Expressing M n P n from Eq. (H3J and denoting ^ — v/v' and 77 = A y< - n ™°)?/, we find 



Using Eq. ( jH]) for Cp™ "*, the right-hand side of Eq. (jl5]l is integrated as 



7* 1 /■ ^OO 

^ A -y P (n-n ) / / ^y Mno P no ( V M0 di drj. 

n =0 JJo 



,(«) 



n-l 
n =0 



yp(n-n ) Jn ) 



where 



'0 



Using Eq. ( 146]) for c p n ^ and c p n 1,1 , we compute the following difference 



and obtain 



c W = A - J /p(i + ^) c C«-D. 



(42) 



^ 1 /*oo 1 

M n P» = V / \y^- n ^M no P no {\y^- n ^v')- ^ (^) dv'. (43) 

n =0 J0 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



Expression (j4"9l implies the power-law (l3"Tj) with the exponents 

C P = l + (p-l)y-log A (l + ^_ 1 ). (50) 

Since <p p -i > 0, this expression satisfies the inequality ( }2"g|) . 

We determined the scaling exponents £ p explicitly in terms of the moments of the creation 
rate function (p(v). For a particular example, when the creation rate is given by self-similarity 



arguments, this computation was done in 



25] . As we showed in Section IVIj if ( p exist and 
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are differentiable for all real values of p, then the Gartner-Ellis theorem ensures the large 
deviation principle fl38|) . Therefore, the large deviation principle in our model is verified for 
any creation rate function, which has finite and differentiable moments (p p for all real p. 

We proved that the large deviation principle emerges naturally for a class of simple models 
of instanton creation. These models have several essential simplifications. In particular, 
we disregarded correlations in time and correlations between instantons in different shells. 
Numerical simulations show that such correlations are important. In fact, numerical values 
for the moments tp p deviate strongly from the values determined by Eq. (1501) for known £ p . 
Also, the function <p(v) depends on a way it is computed. On the other hand, the numerical 
data provided good evidence for universality of the creation process, because the function 
(p(v) found numerically does not depend on the shell number n. 

VIII. TURBULENT REGIME DOMINATED BY A SINGLE BLOWUP 

We showed in Section [V] that the exponent y describing the universal scaling of instantons 
in Eq. (1271) does not determine any of the scaling exponents ( p , neither their asymptotic 
behavior for large p, see Eq. (|28|) and Fig. HI The anomalous scaling of velocity moments 
is linked exclusively to the process of instanton creation. In this section we demonstrate a 
different turbulent regime, where the instanton scaling has strong influence on the velocity 
moments. This regime can be predicted by looking at the blowup scaling exponent yo 
depending on the model parameter in Fig. [2^. At c = —0.139, we have yo = 0. Hence, in the 
asymptotic form of blowup given by Eq. ( 1141) . all the local maxima v n = max t |it n (£)| are of 
the same order of magnitude. As a result, a single blowup provides the terms (Tk n )~ 1 v^ l ~' 1 ~ 
k~ x in the sum f[T9~j) for the structure functions S' , which yields the condition ( p < 1 for 
all p. This upper bound is exact for £3 = 1, which is required by the existence of energy 
cascade. This fact suggests that the blowup scaling should play essential role for the Sabra 
models with the parameter c in the neighborhood of —0.139. 

Figure [TD] (thin lines and circles) shows the structure functions S p (k n ) and S' p {k n ) for 
c = —0.2 in logarithmic coordinates. Vertical shifts are used to compare the graphs for 
S p (k n ) and S' (k n ), and the good match confirms validity of the description based on velocity 
maxima. The inertial range shrinks substantially in this model and corresponds roughly to 
the shells 17 < n < 30. In this case the blowup scaling exponent yo = 0.0534 is small. As 
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FIG. 10. (Color online) Thin black curves present the velocity moments S p (k n ) = (\u n \ p ) for 
p = 1, 3, 5, 7, 9, 11. Circles determine the functions S'(k n ) from Eq. f)19|) for the same p and even n 
(vertical shifts are applied to facilitate comparison with S p ). The red (bold gray) curves show the 
values (Tk n )~ 1 Vn~ 1 for a single dominant instanton and p = 5, 7,9, 11. 

we just mentioned, due to a slow decay of blowup amplitudes in Eq. ( fT4l) . we expect that 
the inertial range is influenced by the blowup scaling. 

The numerical simulation shows that a single stable instanton dominates the inertial 
range. This instanton is created in the initial shell n = and travels all the way to the 
viscous range. The bold red (gray) curves in Fig. [10] show the values of a specific term 
(Tk n )~ 1 v^~ 1 in the sum (119jl . which corresponds to this instanton. In the figure, we used the 
same vertical shift for the red curve as for the full sum S' p (k n ), which shows that not only 
the slope but also the value of S'(k n ) is determined by a single instanton for large p. Due 
to its dominant role, this instanton is weakly influenced by surrounding fluctuations, i.e., by 
other instantons. As a result, we can expect that the dominant instanton scales with the 
same exponent yo as the blowup. This hypothesis agrees perfectly with the numerical data. 

Figure [TT] shows the scaling exponents ( p for the Sabra model with c = —0.2 and —0.139. 
In this case the blowup scaling exponent is equal to yo = 0.0534 and 0, respectively, see also 
Fig. [2^. The dotted straight lines in Fig. [TT] show the right-hand side of the inequality (I28j) 
with y = y . We see that this inequality becomes exact for large p because of the dominant 
role of a single instanton. Note that the graph of ( p is not concave, as it must be, and 

22 



1.6 
1.4 
1.2 
1 



0.6 
0.4 
0.2 







c = -0.2 








c = -0.139 



4 6 8 10 

P 



FIG. 11. (Color online) Anomalous scaling exponent Q p computed for the Sabra model with c = 
—0.2 and —0.139. For each case, the red dotted line shows the values of Cp determined by the 
universal scaling of a single blowup. 

violates slightly Eq. (128]) in the region near the intersection with the dotted line. This seems 
to be a numerical artifact due to very slow convergence in the region where the blowup 
scaling competes with the scaling of instanton creation process. The horizontal part of the 
graph with ( p = 1 for c = —0.139 reminds the analogous behavior of scaling exponents for 
turbulence of the Burgers equation, see 3J, lZ[. Our results show that there are parameter 
values of the Sabra model, when the anomalous scaling exponents ( p are explicitly related 
to the blowup scaling exponent y for some range of p. 



IX. CONCLUSION 

In this paper we have shown that the blowup (a singularity developing in finite time) may 
be considered as a basic element in the theory of developed turbulence for the Sabra shell 
model. We utilize the fact that the blowup in the inviscid system has universal asymptotic 
form, where shell speeds pass successively through their maxima with increasing time and 
wavenumber. This sequence of maxima is the main property used for identifying and analyz- 
ing the coherent turbulent bursts, which are induced by blowup and called the instantons. 

Blowup is characterized by the asymptotically stable traveling wave solution for the 
renormalized system. Thus, almost any initial condition of finite norm leads to blowup in 
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the inviscid Sabra model. In the inertial range of turbulent regime, the blowup-like structures 
appear in every available part of space-time and propagate toward the viscous range. This 
dynamical behavior can be viewed as a "gas" of instantons, which is dense in space-time and 
moves from the forcing to the viscous range. Existence of many interacting instantons alter 
their properties, as compared to the "pure" blowup, but the instanton statistics remains 
self-similar and universal. 

We showed that the anomalous scaling of velocity moments is a natural consequence of 
the instanton creation process, which obeys the large deviation principle. This allows, in 
particular, justifying the universal form of instanton probability density functions in inertial 
range and describing these universal functions analytically in terms of the scaling exponents 
of velocity moments. The obtained results are in excellent agreement with the numerical 
data. 



The described dynamical picture brings us back to the famous Richardson description [301 ] : 
"Big whirls have little whirls that feed on their velocity, and little whirls have lesser whirls 
and so on to viscosity", which is known to be inadequate for the Navier-Stokes turbulence {2]. 
We see now that, in the Sabra shell model, this description becomes true if one substitutes 
"whirls" by "instantons" (or "blowups"). The dissipation anomaly follows naturally from 
this picture, because the instanton dynamics is unrelated to viscosity in the inertial range and 
the instantons move only toward large wavenumbers. The intermittency becomes a simple 
consequence of the instanton scaling, where the scales of time and velocity are related as 
At n ~ (k n v n )~ l . This means that large- amplitude instantons are fast, while small-amplitude 
instantons are slow and can be viewed as windows of low activity. We showed, however, 
that the described scenario is not the only possibility, and the anomalous exponents may be 
linked to universal properties of a single blowup for specific values of the model parameter. 

The essential part of our work is based on numerical data. Here the blowup, whose prop- 
erties follow from the model equations, is used as a guideline for the numerically accessible 
definition of instantons. If the analytical theory of turbulence for the shell model can be 
developed in a similar manner, the formal definition of instantons and derivation of their 
universal scaling properties directly from the shell model equations would be the major dif- 
ficulty. An important step in this direction was done in [26[, where self-similar statistics of 
instantons was derived as a result of the blowup interacting with small random fluctuations. 
Our results show that this theory should be extended by taking into account interactions 
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among instantons and instanton creation. Note that, for different shell models, the blowup 



structure may be not self-similar but still universal 



15] providing an extra complication. It 



would be interesting to test these ideas also on the cascade models of turbulence, where the 
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2l|. 



interactions among variables within the same shell are introduced 

As for the 3D Navier-Stokes developed turbulence, our results confirm the common un- 
derstanding of the importance of coherent structures like, e.g., formation of vortex filaments. 
The novel idea is that it is the universal creation process of these structures what may drive 
turbulent dynamics in the inertial range, while the scaling of an individual structure plays 
a secondary role. Moreover, one may notice that such structures do not have to blow up 
in finite time in the inviscid limit and, e.g., the exponential rate would be sufficient. The 
method for identifying and tracking coherent structures from the moment of their creation 
till the viscous range has to be developed in order to verify our hypotheses numerically or 
experimentally. 
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